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', Abstract 

ON 

It has been recently shown that, in the first order (Palatini) for- 
q-( malism, there is universality of Einstein equations and Komar energy- 

^ 1 momentum complex, in the sense that for a generic nonlinear La- 

grangian depending only on the scalar curvature of a metric and a 
torsionless connection one always gets Einstein equations and Komar's 
expression for the energy-momentum complex. In this paper a similar 
■ analysis (also in the framework of the first order formalism) is per- 

formed for all nonlinear Lagrangians depending on the (symmetrized) 
Ricci square invariant. The main result is that the universality of Ein- 
stein equations and Komar energy-momentum complex also extends 
to this case (modulo a conformal transformation of the metric). 
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1 Introduction 



It has been recently shown |2| that Einstein equations are "univer- 
sal", in the sense that for a generic nonlinear Lagrangian f(R) depending 
only on the scalar curvature R of a metric and a (torsionless) connection 
one always gets Einstein equations if one relies on the first order (" Pala- 
tini ") formalism. The importance of considering nonlinear gravitational 
Lagrangians is particularly stressed by the fact that they provide a simple 
but general approach to governing topology ||; applications to string theory 
have also been considered [fl; for a review see []|, [|. Locally, in fact, one 
gets standard Einstein equations from a nonlinear Lagrangian, but there are 
global topological effects which are described by nonlinear Lagrangians and 
which are absent in the linear case. It is therefore important to investigate 
the universality property not only at the level of the equations of motion, but 
also for the energy of the gravitational field. This problem was addressed by 
us in a previous paper Q , where it was shown that universality holds also for 
the energy-momentum complex. We calculated in the first order formalism 
the covariant energy-density flow for a nonlinear Lagrangian f(R) and we 
showed that, in a generic case, it equals in fact the Komar expression already 
known in the purely metric formalism for the linear Hilbert Lagrangian. 

In this paper we shall extend our discussion to the case of Lagrangians 
with an arbitrary dependence on the square of the symmetric Ricci tensor 
constructed out of a metric and a (torsionless) connection. We shall first 
suitably extend the universality theorems of p], |j concerning field equations, 
by showing that also in this case field equations generically reduce to Einstein 
equations (modulo a conformal change of metric); also the non-generic cases 
will be shortly considered (as in f|, [|]). Subsequently we shall calculate the 
energy-momentum complex for this family of Lagrangians and we shall show 
that it generically reduces to the standard Komar expression, provided the 
metric undergoes a conformal rescaling. In view of further investigations, 
we shall finally provide some calculations concerning functions of the square 
of the full Ricci tensor. 

Our results can be relevant also for quantum gravity. As is well known, in 
fact, to remove divergences in quantum gravity one has to add counterterms 
to the Lagrangian which depend not only on the scalar curvature but also 
on the Ricci and the Riemann tensors invariants. It follows from our results 
that "on-shell", in the first order formalism, counterterms depending on the 
scalar curvature and the Ricci square invariant generically do not change 
the semiclassical limit, since we still have the standard Einstein equations. 
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The main result of the present paper is the following: 

Theorem Let S = g^ g a ^ Ri lxa \(r)Ri v p\(T) be the (symmetrized) Ricci 
square invariant, build out of a metric g and a symmetric connection T. 
For the action 

A{jg,T) = J f(S)Jgdx 

where f is a generic analytic function of one real variable, the Euler-Lagrange 
equations (in the first order formalism) lead always to the Einstein equations 

R^iu(h) = 7/v 

for a new metric h which is defined as h^ u = 7 _1 R^ U ^(T) - with a constant 
7/0. 

Moreover, a superpotential for the energy-momentum flow along any vec- 
torfield £ is given by the Komar expression 

This result will follow from Propositions 2.1, 2.3 and 3.2 below. 

It turns out therefore that the action A(g, T) leads to two metrics: be- 
sides the initial metric g one gets a new Einstein metric h, related to g by 
the algebraic equation 

(g-'hf = ± I 

This condition provides on spacetime some additional structures, namely: 
a Riemannian almost-product structure or an almost-complex structure 
with a Norden metric. These aspects will be studied in more detail in our 
forthcoming paper H. 

2 Field Equations 

As is well known gravitational Lagrangians which are nonlinear in the 
scalar curvature of a metric give rise to higher derivatives or to the appear- 
ance of additional matter fields H This strongly depends on having 
taken a metric as basic variable; moreover the equations ensuing from such 
Lagrangians show an explicit dependence on the Lagrangian itself. It was 
shown in Q that, in contrast, working in the first order (Palatini) formalism 
with independent variations with respect to a metric and a symmetric con- 
nection, then, for a large class of Lagrangians of the form f(R) ^/g, where 
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R is the scalar curvature, the equations obtained are almost independent on 
the Lagrangian, the only dependence being in fact encoded into constants 
(cosmological and Newton's ones). In this sense the equations obtained are 
"universal". These are in fact Einstein equations in generic cases, while 
in degenerate situations and in dimension n = 2 one gets either equations 
which express the constancy of the scalar curvature or conformally invariant 
equations. 

In this paper we shall consider, still in the first order (Palatini) formal- 
ism, the family of actions 

A(T,g)= ( f(S)^gdx (2.1) 
Jm 

where: M is a n-dimensional (connected) manifold endowed with a metric 
g^u of an arbitrary signature and a torsionless (i.e., symmetric) connection 
V a • the Lagrangian density is Lf = f(S)y/g , where / is a given function 
of one real variable, which we assume to be analytic; the scalar S is the 
symmetric part of a Ricci square-invariant, considered as a first order scalar 
concomitant of a metric and the (torsionless) connection, i.e.: 

S = S(g, T) = S^S^ = g^g^SapS^ (2.2) 

where = S^ U {T) = R{nv)(P) is the symmetric part of the Ricci tensor. 
We adopt the standard notation for the Riemann and Ricci tensors of V : 



Rua = Rimt(T) = R^ U(T , (a, ... = 1, ...,n) 

and simply write ^fg instead of | det g^ \ l t 2 . Since the Lagrangian L f does 
not depend on derivatives of the metric field, the only momenta are those 
conjugated to the connection coefficients: 



iwk _ dLf = p dSg/3 



^£j£L = 2(7^ - ir K H?) (2.3) 



where T^ U K = d K T* u and 



i^ = 2y/gf\S)S^ (2.4) 
(recall that = g^ a g vf5 S a p but g^ is the inverse of g^ v ). 
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The Euler-Lagrange equations for the action (2.1) with respect to inde- 
pendent variations of g and T can be written in the following form 

f(S)g^S fia S u(3 - \f{S)g^ = (2.5) 

V a 7r^ = 2VM'(S)y/g Sn = (2.6) 

where V Q is the covariant derivative with respect to T and we assume n > 2. 
In fact, variation of the action (2.1) with respect to T gives the equation: 

V a vr^ - V p ir p ^6$ = 

which, due to the symmetry of tt^ , in any dimension n > 2 reduces to (2.6) 
by taking a trace (see Appendix I). 

The properties of the system of equations (2.5)-(2.6) are governed by the 
following equation for the scalar S: 

77 

f(S)S--f(S) = (2.7) 

which is obtained by transvecting (2.5) with g^ v . Following the lines of the- 
orem 1 of M we shall then distinguish the following three mutually exclusive 
cases: (i) eq. (2.7) has no real solutions; (ii) eq. (2.7) has isolated real so- 
lutions; (hi) eq. (2.7) is identically satisfied. Their discussion proceeds as 
follows: 

(Case 1 — no real solutions) If equation (2.7) has no real solution, then 
also the system (2.5)-(2.6) has no real solution. 

(Case 2 — isolated real solutions) Let us now suppose that eq. (2.7) is 
not identically satisfied and has at least one real solution. In this case, since 
analytic functions can have at most a countable set of zeroes on the real 
line, eq. (2.7) can have no more than a countable set of solutions S = Ci 
(i = 1,2,.. .), where Cj are constants. Consider then any solution 

S = Ci (2.8) 

of eq. (2.7). We have now two possibilities, depending on the value of the 
first derivative /'(cj) at the point S = Ci. 
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(Subcase 2.1 - f'(a) / 0) Let us assume that /'(q) / 0. Then eq. (2.6) 
takes the form 

V a (Vff S^(F)) = (2.9) 
while equation (2.5) reduces to 

g a PS^(T)S u(3 (T) = A( Ci )g^ (2.10) 

where a constant A = A(q) arises according to: 

A = A(cj) = f{ci)/Af{ci) = a/n (2.11) 

We shall discuss separately the subcases n > 2 and n = 2. 

(Subcase 2.1.1 — /'(cj) / and n > 2) As it is known, for n > 2 and any 
metric h^, the general solution of the equation 

Vl(Vdet hh^) = (2.12) 
is the Levi-Civita connection T = r^c(/i) (here is the inverse of /i Mi ,): 

The Ricci tensor i?^(r) is automatically symmetric and in fact identical 
to the Ricci tensor R^ v (h) of the metric h itself. Hence the following holds 
true for the present subcase: 

Proposition 2.1 Assume A / and let us take 7 / and n > 2. Then: 
• (i) If h^ v is a metric satisfying the Einstein equations 

R»u(h) = tV ( 2 - 13 ) 
and a metric g^ is an arbitrary solution of the algebraic equation 

9 aP 9 up K a h v(3 = A$£ (2.14) 

where g^ v denotes the inverse of g^ u , then the pair (g,TLc(h)) satisfies 
eq.s (2.9)-(2.10). 
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• (ii) Conversely, if (g,T) is a solution of eq.s (2.9)-(2.10) and 7 7^ 0, 
then 

has to satisfy the algebraic relations (2.14). Moreover T = Tic{h) and 
one has Einstein equations (2.13). 

Proof. Assume that h and g are two metrics satisfying (2.14). In partic- 
ular this implies that both metrics are nondegenerate and that det h 2 = 
(^) n det g 2 , which ensure proportionality by a constant factor det h ~ det g. 
Moreover one finds 

(h-'r ~ g^g^Kp (A) 

which tells us that rising indices of h by means of g produces a matrix which 
is proportional to the inverse of h. 

Now assume that h is an Einstein metric: Ru, v (h) = jh^ u with 7 / 
and take T = TlcQi). Then 

S^(T) = R, lu (T LC ) = R^(h)= 7 h flu (B) 

Therefore, V ^{^/hh a ^) = implies (2.9). The first claim is thence proved. 

Assume conversely that (2.9) and (2.10) are satisfied. Setting h^ = 
(7 _1 Snv) it then follows from (2.10) that h satisfies the algebraic rela- 
tion (2.14). Now, using (A), equation (2.9) can be rewritten in the form 
\/ a (^/hh fJ,u ) = 0. This in turns implies that T = YlcQi). Therefore, again 
by (A), h is an Einstein metric and our claim is proved. Q.E.D. 



This extends the results previously found by Higgs jjllj] for the confor- 
mally invariant case f(S) = S in dimension n = 4 (see Proposition 2.3), on 
the basis of an earlier paper by Stephenson 12]. The above proposition can 



be re-phrased as follows. If we define a new metric h^ by setting: 

hfiv = 7 

with an arbitrary nonvanishing constant 7, then eq.s (2.9) and (2.10) imply 
that r is the Levi-Civita connection of the new metric h and eq. (2.10) 
reduces to the Einstein eq.s for the metric h 

Rfiv{h) = 7^V 

This, in turn, leads to a constant scalar curvature for the new metric h (in 
fact, it is R(h) = nj) and M will be an Einstein manifold with respect to 
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the new metric h. 

Remark. Notice that for A > eq. (2.14) has an obvious global solution 
of the form g^ u = e~ u) h^ LU with u> = l/2ln(A/^/ 2 ). The case A < is more 
complicated and it will be considered in a forthcoming publication pi. 

(Subcase 2.1.2 — /'(cj) ^ and n = 2) In the case n = 2, instead, 
equation (2.12) allows a further degree of freedom (related with conformal 
invariance) . It has in fact the following general solution ||, |3], [b| : 

= W° u (h,B) = ]^ a {d^h va +d v h^-d a h^) + ^{5ZB^+5lB u -h liy B <T ) 

(2.15) 

where B a is an arbitrary vectorfield and we set B^ = h^ a B a . This is due 
to the fact that only in dimension n = 2 eq. (2.6) cannot be reduced to 
V ' a hav = (since only in two dimensions it does not imply that the Rie- 
mannian volume element of h is covariantly constant along V). A connection 
W(h, B) having the form (2.15) is called a Weyl connection |14]]. Using eq. 

(2.15) , from the definition of S^iT) one has: 

S^(W(h, B)) = ^(R(h) - D a B a )h^ u (2.16) 

where D a denotes the covariant derivative with respect to the metric h^, 
so that eq. (2.13) reduces to the following scalar equation 

R(h) - D a B a = 2 7 (2.17) 

Equation (2.17) is the "universal" equation for 2-dimensional space-times 
PI, |i~3]j; it replaces Einstein equations, which are the "universal equations" in 
dimension n > 2. Equation (2.17) is in fact the equation of constant scalar 
curvature for the metric h and the Weyl connection (2.15), because from 

(2.16) one has: 

Tl(h, B) = h^ u R^{W(h, B)) = R{h) - D a B a (2.18) 

We remark that eq. (2.17) has always infinitely many local solutions, but 
it might have no global analytic solution (depending on the topology of the 
2-dimensional manifold M). In any case, the following holds true for the 
subcase 2.1.2: 
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Proposition 2.2 Assume A ^ 0, 7 ^ and n = 2. Then: 

• (i) If h^v is a metric satisfying equation (2.17) where B a is an arbi- 
trary vectorfield and a metric g^ u is an arbitrary solution of the alge- 
braic eq. (2.14), then the pair (g,W(h, B)) satisfies eq.s (2.9)-(2.10). 

• (ii) Conversely, if (g,T) is a solution of eq.s (2.9)-(2.10), then 

V = 7 _1 5^(r) 

has to satisfy the algebraic relations (2.14). Moreover F = W(h,B) 
for some vectorfield B a and one has eq.s (2.17). 

Proof . The proof goes along the same lines as in Proposition 2.1, with the 
only difference that Einstein eq.s are now replaced by equations (2.17). 

Q.E.D. 

(Subcase 2.2 - f'(a) = 0) Suppose now f'(ci) = 0. Then eq. (2.7) implies 
that also /(cj) = 0, i.e., c% is a zero of order at least two of f(S). In this 
case, eq.s (2.5)-(2.6) are identically satisfied and the only relation between 
g and F is contained in the following equation 

S(g,F)= Ci (2.19) 

This equation represents a genuine dynamical relation between the metric 
and the connection of M, although it is not enough to single out a connection 
r for any given metric g (as it happened, instead, in subcase (2.1) above, 
where F turns out to be the Levi-Civita connection of h and h is algebraically 
related to g). In fact, defining a tensorfield by: 

A^ = A^( 5 ,r) = r^-r^( 5 ) (2.20) 

equation (2.19) can be turned into a quasi-linear first order PDE for the 
unknown A, having the term S(g) — q as a source. The space of solutions of 
this last equation, as functions of the metric g together with its first deriva- 
tives and a number of auxiliary fields, has a complicated structure. In any 
case, it is easy to see that this space contains as a subspace the space of all 
pairs (g, F) satisfying eq. (2.10) for A = Cj/n. We also remark that for Cj = 
this equation covers the case A = which was excluded from Proposition 
2.1 (see above). 
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(Case 3 — the conformally invariant case) We consider now the case in 
which eq. (2.7) is identically satisfied. Under this hypothesis the Lagrangian 
is proportional to: 

f(S) = S n/A (2.21) 

We shall then consider for simplicity the case S > (analogous results will 
be valid for S < and they may extend across S = at least if n = 4A;). 
We stress that for n = 4 this is exactly the linear Lagrangian f(S) = S 
considered in |ll ; p^ ]. 

In this case equations (2.5) and (2.6) read as follows: 

[g aP S^S vP - -S 9lu ,) = (2.22) 
n 

V a (S*?y/gSn = (2-23) 
Notice first of all that, under conformal transformations 

9v» = e^g^ , f = r (2.24) 
where u is an arbitrary function on manifold M, one has 

S = e~ 2uJ S 



Therefore, the action 



A(g,T) = [ S n '^dx 
Jm 



as well as the equations (2.22) and (2.23) are invariant under the transfor- 
mation (2.24), i.e., A(g,f) = A{g,T). 

Also in this case we shall consider separately the two cases n > 2 and 
n = 2. 



(Subcase 3.1 — conformally invariant for n > 2) If n > 2 there are two 
possibilities. If S = everywhere on our manifold then we have only the 
equation 

S(g,T)=0 
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and the discussion proceeds as for eq. (2.19) above. When 5^0 there is 
instead an additional conformal degree of freedom, as noticed earlier in |Tl]] 
and |l2j (and later exploited in |l5j ) . In what follows we restrict ourselves to 
the case when S is a nowhere vanishing function. Then, in fact, eq.s (2.22) 
and (2.23) reduce to the following: 

9 aP S^ a S u p - ^S 9tlu = (2.25) 

V Q (S-ydetS a/3 5n = (2.26) 

In order to get (2.26) we have used (2.25) to calculate first the determinant 
of g as a function of det (Sap) and next replace it into (2.23). The following 
proposition is true: 

Proposition 2.3 Assume S ^ 0, 7 / and n > 2. Then: 

• (i) If h^v is a metric satisfying the Einsten equations 

RfMu(h) = 7/1^ 

where a metric g^ is an arbitrary solution of the algebraic equation 

9 aP g vp h m K p = (2.27) 

with an arbitrary nowhere vanishing function s, then the pair (g,TLc(h)) 
satisfies eq.s (2. 25) -(2. 26). Moreover one has S(g,Fi,c(h)) = s. 

• (ii) Conversely, if (g,T) is a solution of eq.s (2.25)-(2.26), then 

h^u = 7 _1 SW( r ) 

has to satisfy the algebraic relations (2.27) with s = S(g,T). Moreover 
one has T = Tic{h) and R^iji) = 7/1^. 

Proof. Since S is a nowhere vanishing function on a (connected) manifold, 
it is strictly positive or strictly negative. Therefore, it can be written in the 
form S = i e 2uJ . By the conformal change of metric g 1— > e^g, equation (2.25) 
transforms into (2.14) and equation (2.26) transforms into (2.12). (Notice 
that after the conformal transformation S transforms to 1). Then the same 
reasoning as in Proposition 2.1 proves our claim. Q.E.D. 
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Remark. Notice that for s > eq. (2.27) has an obvious global solution of 
the form g^ u = e'^h^ with a conformal factor 2w = ln(As / '717 2 ) . 

We can also restate the result as follows: if (g, T) is a solution of eq.s 
(2.25) and (2.26), then there exists a scalar field ip such that the conformally 
related metric ipg^v satisfies Einstein equations and the connection T is the 
Levi-Civita connection of ipg^v', moreover, the scalar curvature of the origi- 
nal metric g and the connection T equals tp. The origin of this extra scalar 
field ip = e~ u will be discussed elsewhere [ pi in the framework of Legendre 
transformation for metric-affine theories. 



(Subcase 3.2 — conformally invariant and n = 2) If n = 2 equations 
(2.18) and (2.19) simplify to 

9 aP S^ a S v p{Y) - ±S(g, T)g^ = (2.28) 

V a (5-3^5H = (2.29) 
Then the following holds true: 

Proposition 2.4 Assume S 7^ ; 7 7^ and n = 2. Then: 

• (i) If h^ u is a metric satisfying equation (2.17) where B a is an arbi- 
trary vectorfield and a metric gn V is an arbitrary solution of the al- 
gebraic eq.(2.27), with an arbitrary nowhere vanishing function s £ 
F(M), then the pair (g,W(h,B)) satisfies eq.s (2. 28) -(2. 29). More- 
over S(g,W(h, B)) = s. 

• (ii) Conversely, if (g,T) is a solution of eq.s (2.28)-(2.29), then 

has to satisfy the algebraic relations (2.27) with s = S(g,T). Moreover 
r = W(h,B) for some vectorfield B a and one has (2.29). 

Proof. The proof is an obvious combination of those already given for Propo- 
sitions 2.2 and 2.3 . Q.E.D. 
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3 Energy— Momentum Complex 

We are now ready to apply general formulae for the energy-momentum 
complex to the Lagrangian density 

L = L(g^,T^, r^J = f(S) (3.1) 

These general formulae are considered in Appendix II and follow closely the 
earlier works of ours 1 17 , 15 , IS ] ; see also |2^, 21 ] . Our results will extend 
those of our previous paper |?J. 

Our field variables are now g and V, while the only derivatives entering 
the Lagrangian are the first order derivatives of T. Therefore, from (II- 1) 
one has the following definition of the energy-density flow 

E a =p^ a £^ u -L^ (3.2) 

where p^ a are given by (2.3)-(2.4). By using the decomposition (see II— 6) : 

and taking the symmetry of p£ u ^ into account, one has 

where is the covariant derivative with respect to the connection T. From 
this (see (II— 8)) one can easily recognize that: 

T AfT = D rpXap _ apX 

ct u > ^ a Pa 

Substituting these into (11-13) and making use of (2.3)-(2.4) one gets finally 

u^fe) = 2^ CT ^v CT r ] + (3.3) 

The second term of (3.3) vanishes on shell by virtue of (2.6). Therefore we 
have proved the following: 

Proposition 3.1 The energy- density flow (3.2) for the Lagrangian (3.1) 
can be represented on shell under the form = d v U ,M ' , where the superpo- 
tential U^ u is given by: 

Xjv» = 2 Tf]? Vct £"] (3.4) 

Q.E.D. 
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Let us now discuss the expression (3.4) we found for the superpotential. 
It resembles very much Komar superpotential, which is known to be a good 
superpotential in the purely metric formalism |22[ | (Komar superpotential 
is in fact equal also to the Schrodinger-Moller-Mizkievich superpotential; 



for a discussion see, e.g., 23, ^4|]). Our superpotential differs in fact from 
Komar's one by a factor. If the dimension n of space-time is larger than 2 
and we are in a generic position (see Section 2), i.e. if equation (2.4) has a 
solution S = c such that c / and f'(c) / 0, then we know that T = Tic{h) 
for some h given by eq. (2.14). Accordingly, from (2.4) and (2.14) one has 

2 

vr^= 2 1 {\)^y/hh^ (3.5) 



A 



and (3.4) reduces to: 



2 

U" v = 47(21)5 Vfc V^r ] (3-6) 

which is proportional to the Komar expression for h. This corresponds to 
the subcase 2.1. 

At bifurcation points S = c, i.e., in the subcase (2.2), one has instead 
/'(c) = 0. Therefore in this case one has tt^ = and the energy-density 
flow as given by (3.4) is identically vanishing. 

For the conformally invariant Lagrangian L = S n / A ^fg (which in dimen- 
sion n = 4 is just the quadratic Lagrangian S yj~g) one has finally: 

2^ = nS(g,T)^ 

and one can use a conformal transformation h = e u g (for S ^ 0) to reduce 
the system to Einstein equations (this is in fact the subcase 3.1). Otherwise, 
one can decompose the connection V as in (2.20): 

where the tensor A = A (a;) is symmetric in its lower indices. The superpo- 
tential takes then the form 

= 2 vr^ V CT r ] + 2 TT a[fl A±C (3.7) 

where V is the covariant derivative with respect to the metric g. As we 
mentioned above, the tensor A can be eliminated by a suitable conformal 
transformation, so that the relevant part of the superpotential turns out to 
be again proportional, modulo a conformal rescaling of the metric, to the 
Komar superpotential of h. Therefore one gets 
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Proposition 3.2 In the cases described by Propositions 2.1 and 2.3 above 
the superpotential turns out to be proportional to the Komar superpotential 
of h, i.e. 

Q.E.D. 

The two-dimensional case is particularly important (see || [|, [Dill ). I n 
this case T is a Weyl connection with (see (2.7)): 

and the superpotential reduces to: 

~ V^(D^" ] - (3.8) 

This covers subcases 2.1.2 (see Proposition 2.2) and 3.2 (see Proposition 
2.4). 

A discussion of the energy-momentum complex for non-linear Lagrangians 
in the first order formalism, based on the use of a background connection as 



in 1 25], will be discussed elsewhere 



4 Conclusions 

In this paper it was shown that the universality of Einstein equations 
and Komar energy-momentum complex, already found in previous investiga- 
tions about the first order formalism, for nonlinear Lagrangians depending 
on the scalar curvature, extends also to Lagrangians depending on the Ricci 
square invariant. We have preferred to limit ourselves here to discuss in 
detail various important particular cases and emphasize the general picture. 
It seems to us that the universality of Einstein equations and Komar energy- 
momentum complex for nonlinear Lagrangians depending on the scalar cur- 
vature or the Ricci square invariant is not only interesting and important as 
a mathematical results but also in view of furthers applications to concrete 
problems in gravity theories (including cosmological models), both at clas- 
sical and quantum level. Moreover our results find interesting applications 
to signature change in spacetime. We will discuss these matters, as well 
as the algebraic relations between the two metrics g and h in forthcoming 
papers. It will be in particular shown in H that we obtain, in fact, new in- 
teresting geometrical structures on manifold, so called almost-product and 
almost-complex Einsteinian structures. 
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Appendix I 

In this Appendix we apply the formalism to the case in which the scalar 
Q is equal to the full square of Ricci tensor, i.e. when Q = R^ v R^ u . This is 
in fact the case studied by Higgs and Stephenson (ll], Q for the Lagrangian 
L = Q yfg and our result explains the oversimplification made by Higgs 
(following Stephenson). As above, we shall then consider the general class 
of Lagrangians L = Lf = f(Q) y/g, where 

Q(g, T) = g^g^R^{T)R aP {T) (I - 1) 

is the Ricci square-invariant. Instead of (2.3) we have the following expres- 
sion for the momentum 

PT = = = 2 (^A " (I " 2) 

where 

f^ = 2y/gf'(Q)R^ (1-3) 

and the tensor density Tt^ v is no longer symmetric. (Notice that p^ UK main- 
tains instead the symmetry, i.e. pj^ K = 0, due to the symmetry of T). The 
Euler-Lagrange equations with respect to independent variations of g and 
r can thence be written in the following form 

f'(Q)g a(3 (R, ia Ru/3 + R^Rpu) - \f{Q)g^ = (1-4) 

VaTf^ " V p 7T P ^S^ =0 (1-5) 

Accordingly, the "characteristic equation" for (1-4) is given by 

4/'(Q)Q - nf(Q) = (1-6) 
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Now, being tt^ not symmetric, eq. (1-5) cannot in general be replaced by 

VaTf^ =0 (1-7) 



as it was done in [Tl], [L2|]. Of course, any solution of (1-4), (1-7), under 
the additional assumption d a 7t^ a ^ = 0, is a solution of (1-4), (1-5), but the 
general solution is more complicated (see [^7| , |2^] and |^9| for n = 4) . The 



general solution of (1-5) is the following (n > 2): 

= r^(h) + l[S%A p + 6^ - (n - l)A a v) (i 
where /i Q/ g and are defined by 

and 

A° = ^a«7f [/3a] 

(n-l)(n-2)Vfr 

Appendix II 



In this Appendix we recall some formulae related to the energy- 
momentum complex for the first order formalism (see jj], |3(| |3l]] . 

Let us consider on M a set of fields 0*, which we assume to be fields 
of geometric objects over M, together with a Lagrangian density L(0 l ,0* ) 
depending on the first derivatives 0^ = d^(j)\ According to [17, |l^, |lj| the 
canonical energy-density flow E x is defined by the formula 

E x ^E\t)=p ) >C^ l -Le (II -1) 

where £ = £"<9 Q is an arbitrary vectorfield on M, is the Lie derivative 
and the momenta p\ are defined by 

d<f>\ 

The "weak conservation law" 



d\E =0 (II - 3) 



holds along solutions of the equations of motion (here, and in the sequel, 
d\ denotes "total" derivative with respect to x x ). We assume also that the 
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Lagrangian is natural, in the sense that the lift of any diffeomorphism of M 
is a symmetry for L. This means that one has 

C,L ee d a (CL) = —X^ + (II - 4) 

As a consequence, the first variational formula can be written as follows: 

d^ = -^C^ (II -5) 

where 8L/8(f) 1 denotes the variational derivative of L. Expanding the Lie 
derivative C^cp 1 as a linear combination of £ and its covariant derivatives 
with respect to an arbitrary symmetric linear connection on M, i.e. 

erf = & a e + ^v P e + $r v P v CT r (n - 6) 

the first variational formula can be rewritten in the form 

= W (II - 7) 

where 

ee Ei*(S) = + T£ p V p £ a + T% pa V p V a e (II - 8) 



w = w(o = w a e + wy P e + wf?v p v a e (n - 9) 

The coefficients of W vanish on shell because of Euler-Lagrange equations, 
so that (II-3) holds on shell; moreover, they satisfy so called "generalized 
Bianchi identities": 

w P - v,w; + v v v p w? = (ii - 10) 

Notice that in the formulae (II-6) , (II-8) and (11-10), due to the unique- 
ness of the decomposition, one assumes the symmetry conditions for the 
highest coefficients, i.e.: 

$*IH = 0, TgW = 0, WW = 

Therefore, the energy-density flow admits a representation 

E» = & + d p U w (II - 11) 
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where U pp = U pp {^) is a skew-symmetric tensor density (called a "superpo- 
tential") and = E p (^) is the "reduced" energy-momentum flow, which 
vanishes on shell (see The following formulae hold true : 

& = (W£ - V p WDr + W PP V P C (II - 12) 



U pp = (T^d + v v f%M )f* + f£ pu V u C (II - 13) 

where, Tg pv = {A/3)T^ p]v (recall that T% [pu] = 0). 
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